function z = getLm(N,A,alpha,gamma)
aux = gamma*(1-alpha);
z = A.^(aux/(aux+1)).*alpha.^(gamma/(aux+1)).*N.^(1/(aux+1));
end